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Abstract. Within the framework of a simple model of car traffic on a one-lane 
highway, we study the probability for car accidents to occur when drivers do not 
respect the safety distance between cars, and, as a result of the blockage during the 
time T necessary to clear the road, we determine the number of stopped cars as a 
function of car density. We give a simple theory in good agreement with our numerical 
simulations. 



1. Car accidents 

In the past few years, many highway traffic models formulated in terms of cellular 
automata have been studied, both in one [1-5] and two dimensions [6-7] . For the one- 
dimensional case, several topological variations of the basic model have been proposed, 
including road crossing [rU|, road with junction [rT| and two-lane highway [12-14]. 



Recently, experimental features and complex spatiotemporal structures of real traffic 
flows have been investigated [|T^, ^ . 

One of the simplest models is defined on a one-dimensional lattice of L sites with 
periodic boundary conditions. Each site is either occupied by a vehicle, or empty. The 
velocity of each vehicle is an integer between and t'max- If t) denotes the position 
of the ith car at time t, the position of the next car ahead at time t is x{i + 1, t). With 
this notation, the system evolves according to a synchronous rule given by 

x{i,t + l) = x{i,t) + v{i,t+l), (1) 

where 

t>(z, t -|- 1) = min {x{i + 1, t) — x{i, t) — 1, t) — x{i, t — 1) + a, fmax) (2) 

is the velocity of car i at time t + 1. x(z + 1, t) — t) — 1 is the gap (number of empty 
sites) between cars i and i + 1 aX time t, x{i, t) — x{i, t — 1) is the velocity v{i, t) of car 
i at time t, and a is the acceleration, a = 1 corresponds to the deterministic model of 
Nagel and Schreckenberg |jl|, while the case a = fmax has been considered by Fukui and 



Ishibashi |T^. In this last case, the evolution rule can be written 

x{i, t -|- 1) = x{i, t) + min {x{i + l,t) — x{i, t) — 1, t'max)- (3) 

This is a cellular automaton rule whose radius is equal to t'max- The case a < Vmax is 
a second order rule, that is, the state at time t + 1 depends upon the states at times t 
and t — 1. 

We studied the probabihty for a car accident to occur when drivers do not respect 
the safety distance. More precisely, if at time t, the velocity v{i + l,t) of car i + 1 was 
positive, expecting this velocity to remain positive at time t + 1, the driver of car i 
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increases the safety velocity v{i,t + 1) given by by one unit, with a probabihty p. 
The evolution rule (|l]) is then replaced by 

if v{i + l,t)>0, then x{i,t + 1) = x{i,t) + v{i,t + 1) + AV, (4) 

where AV is a Bernoulli random variable which takes the value 1 with probability p 
and zero with probability 1—p. If v{i + 1, t + 1) = 0, it is clear that this careless driving 
will result in an accident. 

When the car density p is less than the critical car density pc = (1 + ^^max)~^, the 
average number of empty sites between two consecutive cars is larger than Vmax, the 
fraction no of stopped cars is zero and no accident can occur. If p > pc, the average 
velocity is less than Wmax, "^o increases with p and careless driving will result in a number 
of accidents. This number will, however, go to zero for p = 1, since, in this case, all 
cars are stopped. The probability for a car accident to occur should, therefore, reach a 
maximum for density p between pc and 1. 

Neglecting time correlations, we may determine an approximate probability for an 
accident to occur. Let n be the number of empty sites between cars i and z + 1 at time 
t. If the three conditions 

0<n<t;„,ax, v{i + l,t)>0, v{i + l,t+l)=0 (5) 

are satisfied then car i will cause an accident at time t+1, with a probability p. Therefore, 
the value Pas of the probability per site and per time step for an accident to occur is 
given by 

Pas = pno(l - no) J2 - pT = W(l - (1 - py-''^')noil - no). (6) 

n=0 

Dividing by p, one obtains the probability per car and per time step for an accident to 
occur 

Pac = pp{l - (1 - p)'^— +^)no(l - no). (7) 

The simplest approximate expression for the fraction no of stopped cars as a function 
of car density, satisfying the above-mentioned conditions, is 

no = . (8) 

1 - Pc 

This approximation is rather crude. In particular, it neglects the fact that no should 
depend upon fmax- However, as shown in Figure |l|, for v^^x = 3 and a = 1, this linear 
approximation is in rather good agreement with our numerical results. Substituting (|^) 
in (1^) yields 

Pa. = pp{l - (1 - pY--^) ^^ ~^'}\'/\ (9) 

I J- Pc) 

Figure |^ represents the probability Pac as a function of p determined numerically and 
its value given by (pi). 
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2. Stopped cars due to blockage 

As a result of an accident (or any other cause such as road works) , the traffic is blocked 
during the time T necessary to clear the road. For a given T, the number N[p) of 
blocked cars is clearly an increasing function of the car density p. To determine the 
expression of A^, we shall distinguish two regimes. 

If P < Pc, the average car velocity is t'max- Since the average number of empty sites 
between two consecutive cars is equal to 

d{p) = ^i^, (10) 

the line of stopped cars increases by one unit during the time interval 

d{p) ^ ^ 

Hence, during the time T the number of blocked cars is given by 

For p = \ the average number of empty sites between two consecutive cars is equal 
to (i(^) = 1, and the average number of blocked cars increases by one unit at each time 
step. N[p) increases, therefore, from T to T + 1 when p increases from p^ to | (note 
that the ffist blocked site is occupied). When p > \, d{p) < 1, and the average number 
of blocked cars is then given by 

Nip) = T+-l- = T+-^. (13) 
d{p) 1 - p 

Figure ^ represents the average number of blocked function of the car density 

p. The agreement with the approximate expressions of N{p), given by (12) and (0), is 
very good. 
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Figure 1. Fraction of stopped cars hq as a function of the car density p for Umax — 3 
and a — I. The solid hne represents the hnear approximation uq = {p — pc) / {1 — Pc)- 




Figure 2. Probabihty Pac per car and per time step for an accident to occur as a 
function of the car density p for timax — 3 and a — 1. The sohd hne corresponds to the 
approximation given by (H). 



7 



N{p) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P 



Figure 3. Average number of blocked cars as a function of the car density p. The 
sohd hne corresponds to the approximate expressions dl2) and (13). 



